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Abstract 

We measure the subjet multiplicity M in jets reconstructed with a successive 
combination type of jet algorithm (k±). We select jets with 55 < pr < 100 
GeV and \n\ < 0.5. We compare similar samples of jets at yfs = 1800 and 
630 GeV. The HERWIG Monte Carlo simulation predicts that 59% of the 
jets are gluon jets at ^fs = 1800 GeV, and 33% at = 630 GeV. Using this 
information, we extract the subjet multiplicity in quark (M q ) and gluon (M g ) 
jets. We also measure the ratio R = |^ g |_j = 1-84 ± 0.15(stat)l^g(sys). 
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I. INTRODUCTION 



The Tevatron proton-antiproton collider is a rich environment for studying high energy 
physics. The dominant process is jet production, described in Quantum Chromodynamics 
(QCD) by scattering of the elementary quark and gluon constituents of the incoming hadron 
beams. In leading order (LO) QCD, there are two partons in the initial and final states 
of the elementary process. A jet is associated with the energy and momentum of each 
final state parton. Experimentally, however, a jet is a cluster of energy in the calorimeter. 
Understanding jet structure is the motivation for the present analysis. QCD predicts that 
gluons radiate more than quarks. Asymptotically, the ratio of objects within gluon jets to 
quark jets is expected to be in the ratio of their color charges Ca/Cf = 9/4 [|TJ. 

II. THE k T JET ALGORITHM 

We define jets in the D0 detector with the k± algorithm 0-§J. The jet algorithm 
starts with a list of energy preclusters, formed from calorimeter cells or from particles in a 
Monte Carlo event generator. The preclusters are separated by ATI = y/ Ai] 2 + Acf) 2 > 0.2, 
where r\ and are the pseudorapidity and azimuthal angle of the preclusters. The steps of 
the jet algorithm are: 

1. For each object i in the list, define da = where px is the energy transverse to 

the beam. For each pair (i, j) of objects, also define = ^n^^p^j)-^, where D is a 
parameter of the jet algorithm. 

2. If the minimum of all possible da and d^ is a d^, then replace objects i and j by their 
4-vector sum and go to step 1. Else, the minimum is a da so remove object i from the list 
and define it to be a jet. 

3. If any objects are left in the list, go to step 1. 

The algorithm produces a list of jets, each separated by ATI > D. For this analysis, 
D = 0.5. 

The subjet multiplicity is a natural observable of a jet [[|-|9|,|4],§] . Subjets are defined 
by rerunning the kj_ algorithm starting with a list of preclusters in a jet. Pairs of objects 
with the smallest dy are merged successively until all remaining > ycutPr{j€-t). The 
resolved objects are called subjets, and the number of subjets within the jet is the subjet 
multiplicity M. The analysis in this article uses a single resolution parameter y cut = 10~ 3 . 

III. JET SELECTION 

In LO QCD, the fraction of final state jets which are gluons decreases with x ~ px/ ' \fs, 
the momentum fraction of initial state partons within the proton. For fixed px, the gluon 
jet fraction decreases when yfs is decreased from 1800 GeV to 630 GeV. We define gluon 
and quark enriched jet samples with identical cuts in events at ^fs = 1800 and 630 GeV to 
reduce experimental biases and systematic effects. Of the two highest pr jets in the event, 
we select jets with 55 < pr < 100 GeV and \r]\ < 0.5. 
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IV. QUARK AND GLUON SUBJET MULTIPLICITY 



There is a simple method to extract a measurement of quark and gluon jets on a statistical 
basis, using the tools described in the previous sections. M is the subjet multiplicity in a 
mixed sample of quark and gluon jets. It may be written as a linear combination of subjet 
multiplicity in gluon and quark jets: 

M = fM g + (l-f)M q (4.1) 

The coefficients are the fractions of gluon and quark jets in the sample, / and (1 — /), 
respectively. Consider Eq. flOp for two similar samples of jets at yfs = 1800 and 630 GeV, 
assuming M g and M q are independent of \fs. The solutions are 

yl800^j630 _ J630j^l800 
M 1 = J1800 _ J630 ( 4 - 2 ) 

(1 - / 630 ) M 1800 - (1 - / 1800 ) M 630 , . 

9 = 1 ) /180 o_>30 ) ( 4 - 3 ) 

where M 1800 and M 630 are the experimental measurements in the mixed jet samples at 
y/s = 1800 and 630 GeV, and / 1800 and f 630 are the gluon jet fractions in the two samples. 
The method relies on knowledge of the two gluon jet fractions. 



V. RESULTS 
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Subjet multiplicity M 

FIG. 1. Raw subjet multiplicity in fully simulated Monte Carlo quark and gluon jets, 
visibility, we shift the open symbols horizontally. 



For 



The HERWIG 5.9 [|IIJ Monte Carlo event generator provides an estimate of the gluon jet 
fractions. The method is tested using the detector simulation and CTEQ4M PDF. We tag 
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every selected jet in the detector as either quark or gluon by the identity of the nearer (in 
7] x space) final state parton in the QCD 2-to-2 hard scatter. Fig. [I] shows that gluon jets in 
the detector simulation have more subjets than quark jets. The tagged subjet multiplicity 
distributions are similar at the two center of mass energies, verifying the assumptions in 

We count tagged gluon jets and find / 1800 = 0.59 ± 0.02 and / 630 = 0.33 ± 0.03, where 
the uncertainties are estimated from different gluon PDF's. The nominal gluon jet fractions 
and the Monte Carlo measurements at y/s 



1800 and 630 GeV are used in Eqs. (p|~^3|). 
The extracted quark and gluon jet distributions in Fig. [I] agree with the tagged distributions 
and demonstrate closure of the method. 
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FIG. 2. Raw subjet multiplicity in jets from D0 data at -^/s = 1800 and 630 GeV. 



Figure |2] shows the raw subjet multiplicity in D0 data at t/s = 1800 GeV is higher than 
at a/s = 630 GeV. This is consistent with the prediction that there are more gluon jets at 
■s/s = 1800 GeV compared to yfs = 630 GeV, and gluons radiate more than quarks. The 
combination of the distributions in Fig. ||] and the gluon jet fractions gives the raw subjet 
multiplicity distributions in quark and gluon jets, according to Eqs. Ql.2H4.3j ). 

The quark and gluon raw subjet multiplicity distributions need separate corrections for 
various detector-dependent effects. These are derived from Monte Carlo, which describes 
the raw D0 data well. Each Monte Carlo jet in the detector simulation is matched (within 
ATI < 0.5) to a jet reconstructed from particles without the detector simulation. We tag 
detector jets as either quark or gluon, and study the subjet multiplicity in particle jets M ptcl 
vs. that in detector jets M det . The correction unsmears M det to give M ptcl , in bins of M det . 
Figure [3] shows the corrected subjet multiplicity is clearly larger for gluon jets compared to 
quark jets. 

The gluon jet fractions are the largest source of systematic error. We vary the gluon 
jet fractions by the uncertainties in an anti-correlated fashion at the two values of \fs to 
measure the effect on R. The systematic errors listed in Table | are added in quadrature to 
obtain the total uncertainty in the corrected ratio R = f^/lzl — 1-84 ± 0.15(stat)^g^g(sys). 
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FIG. 3. Corrected subjet multiplicity in quark and gluon jets, extracted from D0 data. 



Source 


SR 


Gluon Jet Fraction 


+0.17 
-0.10 


Jet pt cut 


±0.12 


Detector Simulation 


±0.02 


Unsmearing 


±0.07 



TABLE I. Systematic Errors. 



VI. CONCLUSION 

We extract the y cu t = 10~ 3 subjet multiplicity in quark and gluon jets from measurements 
of mixed jet samples at y^s = 1800 and 630 GeV. On a statistical level, gluon jets have more 
subjets than quark jets. We measure the ratio of additional subjets in gluon jets to quark 
jets R = 1.84 ± 0.15(stat)^Q^g(sys). The ratio is well described by the HERWIG parton 
shower Monte Carlo, and is only slightly smaller than the naive QCD prediction 9/4. 
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